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ABSTRACT 

Context. In the last few years we have developed stellar model atmospheres which included effects of anomalous abundances and 
strong magnetic field. In particular, the full treatment of anomalous Zeeman splitting and polarized radiative transfer were introduced 
in the model atmosphere calculations for the first time. The influence of the magnetic field on the model atmosphere structure and 
various observables were investigated for stars of different fundamental parameters and metallicities. However, these studies were 
purely theoretical and did not attempt to model real objects. 

Aims. In this investigation we present results of modelling the atmosphere of one of the most extreme magnetic chemically peculiar 
stars, HD 137509. This Bp SiCrFe star has the mean surface magnetic field modulus of about 29 kG. Such a strong field, as well as 
clearly observed abundance peculiarities, make this star an interesting target for application of our newly developed model atmosphere 
code. 

Methods. We use the recent version of the line-by-line opacity sampling stellar model atmosphere code LLmodels, which incorporates 
the full treatment of Zeeman splitting of spectral lines, detailed polarized radiative transfer and arbitrary abundances. We compare 
model predictions with photometric and spectroscopic observations of HD 137509, aiming to reach a self-consistency between the 
abundance pattern derived from high-resolution spectra and abundances used for model atmosphere calculation. 
Results. Based on magnetic model atmospheres, we redetermined abundances and fundamental parameters of HD 137509 using 
spectroscopic and photometric observations. This allowed us to obtain a better agreement between observed and theoretical parameters 
compared to non-magnetic models with individual or scaled-solar abundances. 

Conclusions. We confirm that the magnetic field effects lead to noticeable changes in the model atmosphere structure and should be 
taken into account in the stellar parameter determination and abundance analysis. 
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1 . Introduction mentioned above and to answer the question how does the mag- 
netic field act at different temperatures and what one could ex- 

The atmospheric structure of magnetic chemically peculiar (CP) pect if the magnet ic field is ignored in calculations of model at- 

stars deviates from that of normal stars with similar funda- mosphere f magnetic CP stars. It was shown that the strength of 

mental parameters due to unusual chemistry, abundance mho- the magnet ic field is the key characteristic controlling the mag- 

mogeneities and the presence of strong magnetic field. These ^tude of the magnetic field effects and the polarized radiative 

effects are not considered in the standard model atmosphere transfer should be taken int0 account . In contrastj me orienta- 

calculations, possibly leading to errors in the stellar parame- tion of the magnet ic field vector does not have much influence 

ter determination and abundance analysis. To circumvent this on any of the obserV ed stellar characteristics and, thus, can be 

long-standing problem of stellar astrophysics, we have devel- safety i gnore( j j n t h e analysis routines, 
oped a new lin e-by-line opacity s ampling model atmosphere 

code LLmodels (Sh ulvak etal.ll200l . Using this tool in the se- So far our models with the magnetic field effects included 

ries of recent papers, we in vestigated in detail the effects of have been developed and applied only in the context of purely 

anomalous Zeema n splitting (IKochukhov et al.l l2005h. polarized theoretical studies. Here we make the first attempt to model the 



radiative transfer dKhan~& Shulyak 2006a) and inclination of the atmosphere of a real star. In this work we use the LLmodels stel- 
magnetic field vector dKhan & Shulva k 2006b) on the model lar model atmosphere code to investigate the atmospheric struc- 



structure, energy distribution, hydrogen line profiles, photomet- ture of the star > takin 8 int0 account individual chemical composi- 

ric colors and the magnitude of bolometric corrections for a grid tlon ' anomalous Zeeman splitting and polarized radiative trans- 

of model atmospheres with different effective temperatures and ^ er - 

metallicities. For the first time we were able to obtain new re- HD 137509 (HIP 76011, NN Aps) is a B9p SiCrFe chemi- 

sults applying direct and self-consistent modeling of all effects cally peculiar star with a strong reversing longit udinal field and 

variable lines of Sin and iron-peak elem ents dMathvsl Il991h 

Send offprint requests to: D. Shulyak, iMathvs & Lairzl Il997l) . iKochukhovl d2006l) (hereafter Paper I) 

e-mail: denis@jan. astro. univie.ac. at has detected resolved Zeeman split lines in the spectrum of 
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HD 137509, showing that this star is characterized by a non- 
dipolar magnetic field geometry with a mean surface field 
strength of about 29 kG. This is the second-largest magnetic field 
ever found in a CP sta r (the first plac e is occupied by the well- 
known Babcock's star. |Prestonl(ll969l) ). The atmospheric param- 
eters, r eff = 12 750 ±50QK and log g = 3.8 + 0.1, were derived 
in Paper I using theoretical fit t o the observed H/3 and Hy line 
profiles based on the ATLAS9 dKuruczlll993al) model with en- 
hanced metallicity, [M/H] = +1.0. The appearance of such a 
strong magnetic field and the presence of outstanding abundance 
anomalies inferred in Paper I allow us to use HD 137509 as a test 
ground for the application of the new generation magnetic model 
atmospheres. 

In the next section we briefly describe the techniques em- 
ployed to construct magnetic model atmospheres. Results of the 
calculations for HD 137509 are presented in Sect. [3] Main con- 
clusions of our study are summarized in Sect. [4] 



2. Model atmosphere calculations 

To calculate magnetic model atmospheres we used the cur- 
rent version of the LLmodels code, originally developed by 
IShulvak et all d2004l) . The LLmodels is a LTE, 1-D, plane- 
parallel, hydrostatic stellar model atmosphere code for early and 
intermediate-type stars. The direct line-by-line calculation of the 
bound-bound opacities implemented in this code allows one to 
account for arbitrary individual and stratified stellar abundance 
patterns and include various effects caused by t he magnetic 
field The VALD database dPiskunov et alj 119951: iKupka et al.l 
1 19991) was used as a source of atomic line parameters. The ex- 
tracted lines were subjected to a preselection procedure inside 
LLmodels, with the default criterion an ne /a cont > 1% in at least 
one atmospheric layer allowing to select only those lines that 
significantly contribute to the opacity (here ai; ne and a cont are 
the line center and continuum opacity coefficients, respectively). 
The numerical experiments showed that this criterion is suffi- 
cient for accurate representation of the energy distribution and T- 
P structure of the models. Decreasing this value leads to signifi- 
cant increase of the number of selected lines but w ithout produc- 
ing n oticeable changes in overall blanketing (see Shulvak et al. 
2004). The preselection procedure was performed twice per 
model atmosphere calculation: at the first iteration and at the iter- 
ation when the temperature correction at each model atmosphere 
layer is less than few tens of K. The remaining model iterations 
are performed with the later line list, thus ensuring a consistency 
of the preselected lines and the model structure obtained. The 
spectrum synthesis for the magnetic model atmosphere calcula- 
tions was carried out in the range between 100 and 40 000 A, 
with a constant wavelength step of 0. 1 A. 

Generally, due to the variation of magnetic field vector across 
the visible stellar surface, one has to compute a number of local 
model atmospheres for each appropriately chosen surface grid 
element using individual values of the strength and inclination 
of the magnetic field. Then, the total flux coming from the star 
should be obtained by integration of the radiation field intensity 
correspon ding to individua l surfac e zones. However, as it was 
shown by Kha n & Shulvakl d2006bl) . the inclination of the mag- 
netic field vector does not significantly influence the structure of 
magnetic models and the resulting energy distribution, i.e. the 
anisotropy effects can be neglected in the magnetic model at- 
mosphere calculations. Consequently, we assume the magnetic 
field vector to be perpendicular to the atmosphere normal and, 
according to the results of Paper I, adopt the field strength of 



Table 1. Elem ent abundances of HD 137509 compared to the 
solar values (Aspl und et al.l 120051) . The second and the third 
columns give, respectively, abundances derived using the scaled- 
solar abundance model atmosphere (Paper I) and using model 
with T e ff = 13750K, \ogg = 4.2 with magnetic field included 
(see Sec. 3.4). Abundances are given in the logarithmic scale 

log(AWNtotal)- 



Element 


tl2750g3.8 (Paper I) 


tl3750g42 


Sun 


He 


< -3.50 


< -3.50 


-1.10 


Si 


-3.73 


-3.58 


-4.53 


Fe 


-3.19 


-3.00 


-4.59 


Cr 


-4.20 


-3.90 


-6.40 


Ti 


-4.54 


-4.20 


-7.14 


Ca 


-7.93 


-7.50 


-5.73 


Mg 


-5.71 


-5.50 


-4.51 



29 kG. The approach for calculating magnetic models used here 
is equ ivalent to the one described in Sect. 2 in lKhan & Shu lvak 
d2006al) . 

For all models presented in this study we assume a force- 
free configuration of the surface magnetic field. This agrees 
with the results o f multipolar modeling of magnetic topology 
dKochukhovl2006l) and implies that possible mod ification of the 
hydros tatic equilibrium by the Lorentz force (see Shulva ket al.l 
d2007l) and references therein) is absent in HD 137509. Thus, all 
differences in the pressure structure between magnetic and non- 
magnetic models that we obtain are caused only by additional 
opacity in the Zeeman components. 

The necessity to perform polarized radiative transfer calcu- 
lations over a wide wavelength range makes the models with 
magnetic field very computationally expensive. To reduce com- 
putational costs, we used the following approach. Fixing funda- 
mental parameters and abundance pattern we calculated a non- 
magnetic model. After this model is converged, a model with 
magnetic field is iterated using the temperature-pressure struc- 
ture of the initial non-magnetic model as a first approximation. 

The abundances of seven chemical elements for HD 137509 
were derived in Paper I and are listed in TableQ]second column). 
Due to the abnormal weakness of the He lines, the value of the 
He abundance should be considered as an upper limit only. One 
can note the strong overabundance of Si, Fe, Cr and Ti, whereas 
Ca and Mg are underabundant by more than 1 dex relative to the 
solar chemical comp osition. We assume solar abundances from 
Aspl und et al.l d2005l) for all other elements. Using available ob- 
servations it is not possible to derive abundances of elements 
other than those listed in Table Q] with good accuracy due to 
complex line shapes and strong blending by iron-peak elements. 
There are no usable lines of intrinsically abundant elements such 
as C, N and Ne, while the abundance of O that could in princi- 
ple be determined from the IR triplet (7772-7775 A) is unreli- 
able due to large NLTE effects and instrumental artifact in the 
UVES spectrum of HD 137509 in this regio n. According to our 
recent investigation dKhan & Shulvakl2007l) . C, N, and O belong 
to the group of elements that produces only a small change in the 
model structure of A and B stars even if their abundances are re- 
duced by 1 dex relative to the solar values. Finally, to ensure 
the consistency between the model atmosphere and the abun- 
dance pattern of the star, we constructed stellar atmosphere mod- 
els with these individual abundances including magnetic field. 
Then we try to assess relative importance of this sophisticated 
approach. 
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Fig. 1. Difference in temperature (upper panel) and gas pressure 
(bottom panel) of the models calculated with individual abun- 
dances (solid line) and with individual abundances +magnetic 
field (dashed line) with respect to the reference model atmo- 
sphere computed with scaled-solar abundances ([M/H] = +1). 
For all models r e ff = 12 750 K and logg = 3.8 is adopted. 



3. Results 

3.1. Model structure 

The effect of the peculiar abundance pattern and magnetic field 
on the temperature-pressure structure of the model atmosphere 
of HD 137509 is shown in Fig. Q] The model with individual 
abundances exhibits a decrease of temperature in the surface lay- 
ers, while the layers close to the photosphere (log T5000 ~ 0) are 
heated compared to the reference scaled-solar abundance model 
([M/H] = +1.0). The effect of magnetic field is, therefore, con- 
siderably different compared to a change in metal abundance. 
The inclusion of magnetic field leads to heating of the surface 
layers due to additional opacity in the Zeeman components. This 
occurs because a high line density in the dominant UV part of the 
stellar spectral energy distribution makes it possible for the back- 
warming effect to occur eve n in surface layers whi ch become 
opaque at these wavelengths (Kochukh ov et alj |2005). The tem- 
perature distribution of the magnetic model in the surface layers 
tends to be close to that of the scaled-solar model. At the same 
time, magnetic field appears to influence the temperature struc- 
ture more efficiently in deep atmospheric layers compared to the 
non-magnetic model with individual abundances. 

Both magnetic and non-magnetic models show a decrease of 
the gas pressure throughout the atmosphere compared to the ref- 
erence model. However, the model with the magnetic field shows 
a steeper depth dependence of the pressure difference relative to 
the model where only individual abundances are taken into ac- 
count. In our hydrostatic models the decrease of the gas pressure 
is caused by the respective increase of the radiative pressure due 



4890 



Fig. 2. Observed and calculated HB line profiles. Upper panel: 
thick line - UVES observations (Paper I), thin line - the log g = 
3.8 model with solar-scaled abundances [M/H] = +1, dashed 
line - the model with individual abundances, dotted line - 
the model with both individual abundances and the magnetic 
field. Bottom panel: difference (in per sent) between model 
with only individual abundances (solid line) and with individual 
abundance+magnetic field (dotted line) relative to the scaled- 
solar abundance model (T e g = 12 750K, \ogg = 3.8 for all 
models). 



to enhanced opacity for models with individual abundances and 
magnetic field. Note that the model with the magnetic field in- 
cluded exhibits a noticeable increase of the radiative pressure in 
the outer atmospheric layers compared to other models. This be- 
haviour of the radiative pressure and corresponding changes in 
the overall temperature-pressure structure of the magnetic atmo- 
sphere are relevant for modern studies of the radiatively-driven 
chemical diffusion in the atmospheres of magnetic stars (e.g., 
lAlecian & Still2007h . 

At this point it is difficult to predict the total effect of the 
magnetic field on the spectral characteristics of the star. In two 
different atmospheric regions (surface and photospheric layers), 
the magnetic field influences the model structure in a differ- 
ent manner, so that the properties of magnetic atmosphere of 
HD 137509 are not equivalent to a usual non-magnetic model 
with a different set of fundamental parameters. One can expect 
that spectral features which are selectively sensitive to either 
temperature or pressure may show different behavior. In the fol- 
lowing section we examine the overall effect of the individual 
abundances and magnetic field on the hydrogen B aimer line pro- 
files and metallic line spectra. 
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Fig. 3. Same as in Fig.[2]but for the Hy line. 



3.2. Hydrogen line profiles and metallic line spectra 

The introduction of peculiar abundances taken from Paper I (see 
Table [T) to the model atmosphere calculations for HD 137509, 
together with the magnetic field, leads to noticeable changes in 
the hydrogen line profiles. In Figs. [2] and [3] we compare theo- 
retical profiles for different model atmospheres with those ob- 
served for HD 137509. Synthetic profiles w ere calculated using 
the Synthmag program ( IKochukhovll2007l) . This code incorpo- 
rates rec ent improvements in the treatment of the hydrogen line 
opacity (Barkl em et al.l l2000) and takes into account magnetic 
splitting of hydrogen lines. However, poss ible modification of 
the Stark broadening by magnetic field (Math vs et al.ll200()f) is 
neglected. The observed profiles of hydrogen lines are extracted 
from the UVES spectrum of HD 137509 described in Paper I. 

One can see that in order to retain acceptable agreement be- 
tween observations and theory, it is necessary to increase the 
logg value from 3.8 derived in Paper I to 4.0. Note, that most 
of the changes in the hydrogen line profiles are due to anoma- 
lous abundances adopted for the star (mainly due to extreme He 
underabundance). Nevertheless, the changes in the temperature- 
pressure structure of the atmosphere produced by the magnetic 
field are also important and has to be incorporated to the model 
in order to retrieve an accurate estimate of the gravitational ac- 
celeration. 

The metallic line spectra are not particularly sensitive to the 
changes in the atmospheric model structure associated with the 
inclusion of the magnetic field. Fig. [4] shows the observed and 
synthetic profiles of some prominent spectral lines of silicon. 
All theoretical spectra were calculated using the Synthmag code 
with the same abundances and the homogeneous surface mag- 
netic field distribution, characterized by (B) = 29 kG, but us- 
ing two different model atmospheres: one model is only with 
individual abundances and another one includes both anomalous 




Fig. 4. Influence of magnetic model structure on the profiles of 
selected spectral lines for T e g = 12 750 K, logg = 4.0 models. 
Double line - UVES observations, thin line - spectra calculated 
using the non-magnetic atmospheric model, dashed line - with 
the magnetic field included. 



abundances and the magnetic field. For both models the funda- 
mental parameters are T e g = 12 750 K, logg = 4.0. The calcu- 
lated spectra were convolved with v sin i = 20 km s _I rotational 
broadening to allow comparison with observational data. The ef- 
fects of the magnetic field in the model atmosphere distort dif- 
ferent lines in different ways. Some of the spectral features con- 
sidered here are not sensitive to the magnetic model atmosphere 
effects at all (for example, Sin 41 30.872 A and Sin 41 30.894 A), 
while other lines show a noticeable discrepancy between the 
spectra computed for magnetic and non-magnetic atmospheres 
(e.g., Tin 4129.161 A). The average difference between the line 
profiles obtained from the magnetic and non-magnetic models 
was found to be about few per sent for the spectra already broad- 
ened by rotation. Due to a very large magnetic broadening (com- 
parable to the rotational Doppler effect), the effect on uncon- 
volved spectra is roughly of the same magnitude for all spec- 
tral lines considered. These results are in a good agre ement with 
our previous investigation dKochukhov et al.l 120051) . However, 
the present results are more accurate because polarized radiative 
transfer was taken into account for both the model atmosphere 
and the spectrum synthesis calculations. 



3.3. Energy distribution 

Generally, for self-consistent modelling of stellar atmospheres, 
the following observables must be reproduced simultaneously: 
hydrogen line profiles, energy distribution and metallic line 
spectra. Among these, the flux distribution is especially sensi- 
tive to the overall energy balance in the stellar atmosphere over 
a wide range of optical depths. So, this observable is extremely 
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useful for determining basic stellar parameters and is prefer- 
able in this respect to broad- and intermediate-band photometric 
data. The role of energy distributions is even more important for 
chemically peculiar stars because most of the widely used pho- 
tometric calibrations are based on observations of normal stars 
and, hence, may not be applica ble to stars with strong magnetic 
fields and unusual abundances dKhan & Shulvak|[2007l) . 

Figure [5] illustrates the effects of individual abundances and 
the magnetic field on the energy distributions for the models 
of HD 137509 with 7/ eff = 12 750K, logg = 3.8. It is evident 
that both the model with individual abundances and the model 
combining the effects of unusual chemical composition and the 
strong magnetic field exhibit flux redistribution from the UV to 
visual wavelength region. The most significant effect is found for 
the model with the magnetic field included, which is in agree- 
ment with our previous results. 

Due to the increased level of flux in the visual region, the 
models with magnetic line blanketing are characterized by an 
anomalous bolometric correction. We find that the difference 
in this parameter with respect to the computations done for 
solar chemical composition reaches ABC = 0.1 mag for the 
model with individual abundance, whereas ABC « 0.15 mag 
for the model with individual abundance and magnetic field. 
Thus, the anomaly in the bolometric correction is quite substan- 
tial and has to be taken into account in the determination of ab- 
solute luminos ity and comparison with evo lutionary models of 
CP stars (e.g., Kochukhov & Bagnulo 2006), yet, even for such 
an extreme star as HD 13750 9, we cann ot confirm the reality of 
ABC > 0.2 mag proposed bv lLanzl (11984 . 

Unfortunately, no suitable observed energy distribution for 
HD 137509 which could be used for verifying theoretical mod- 
els is available in the literature. One can possibly use the flux- 
calibrated spectrum from the UVES LibrarjfJ T he description 
of the reduction of these data can be found in Ba gnulo et al.l 
(2003). Despite the fact that the UVES pipeline provides a user 
with the spectra calibrated in relative units, these data can not 
be used for comparisons with models in a wide spectral range 
due to the flux calibration uncertainties (S. Bagnulo, private 
communication) that could affect the shape of the energy dis- 
tribution. Nevertheless, the relative shape of the UVES spectra 
within small spectral regions (300-500 A) is reasonably well- 
determined. Therefore, in Fig.[6]we compare relative fluxes (nor- 
malized at A — 5000 A) of the magnetic and non-magnetic mod- 
els of HD 137509 with the UVES spectrum in two short regions. 
Note, that both theoretical models presented in the figure were 
calculated with the abundances from Paper I. Thus, the differ- 
ence between the two model fluxes is entirely due to the mag- 
netic field, which is responsible for producing complex spectral 
features clearly seen in the observed spectra. In particular, the 
Zeeman splitting influences the strength of the line absorption 
around A - 5200 A, which is associated with the well-known 
depression in the spectra of CP stars and its photometric char- 
acteristic (Aa photometry, discussed below). Note, that for the 
comparison of fluxes in the UV region (upper panel in Fig. [6]l we 
had to shift the energy distribution produced by non-magnetic 
model along the y-axis to match the observed data. This is a con- 
sequence of the substantial difference in the fluxes close to the 
Balmer jump where the model with the magnetic field shows a 
prominent flux excess compared to the non-magnetic model (see 
Fig.©. 

The only other possibility would be to use the flux cali- 
brated spectra of HD 137509 obtained by the Far Ultraviolet 
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Fig. 6. Comparison between the observed and calculated fluxes 
of HD 137509 in the near UV (upper panel) and visual region 
(lower panel). Gaussian smearing with resolution R = 1500 is 
applied to all the spectra to provide a better view. In the upper 
plot the fluxes for the non-magnetic model are shifted along the 
vertical axis to match the observed spectrum. 

Spectroscopic Explore^ (FUSE) mission. However, these data 
were obtained in a very short wavelength range (910-1 185 A), 
which makes it impossible to use these spectra for accurate de- 
termination of fundamental stellar parameters and testing theo- 
retical models. For example, the overall radiative energy emit- 
ted inside 910-1 185 A region is about 1.6% of the total flux 
for the model with T e ff = 12750K and about 2.5% for the 
?eff = 13750K model, which is unimportant for overall radia- 
tive energy balance in the atmosphere of HD 137509. Moreover, 
the slope of the energy distribution in this region remains the 
same for both models mentioned above. This means that, based 
on FUSE data, it is impossible to distinguish models different by 
as much as 1000 K. Well-calibrated energy distribution covering 
a wide spectral regions is needed. 

3.4. Photometric colors 

As a next step, we have calculated a grid of model atmospheres 
with different effective temperatures (T e s = 12 250K-13 750K, 
AT = 250 K) and gravities (logg = 3.8-4.2, Alogg = 0.2) to 
assess the ability of our models to reproduce the observed photo- 
metric properties of HD 137509. To investigate the role of mag- 
netic field, we have included in the grid both magnetic and non- 
magnetic models with individual abundances. We also consid- 
ered the reference scaled-solar abundances model from Paper I. 

The theoret ical colors were calculat ed following the proce- 
dure outlined in Kochukho v et al.l (120051) . We use modified com- 
puter codes by iKuruczl dl993al) . which take into account trans- 
mission curves of individual photometric filters, mirror reflec- 
tivity and a photomultiplier response function. The synthetic Aa 
values were computed with respect to the theoretical normal- 
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Fig. 5. Theoretical energy distributions for the models with T e ff = 12 750 K, log g = 3.8. Solid line - scaled-solar abundance model, 
dotted line - the model with individual chemical composition, dashed line - the model with individual abundances and the magnetic 
field. Energy distributions are convolved with a Gaussian profile with FWHM = 15 A. 



ity line, ao, determined in Khan&SMy^d2007|). The redden- 
ing, corresponding to E(B - V) = 0.06, was taken into account 
for a ll color indices, except X , Y and Z which are reddening- 
free dc ramer &Maederi[l979l) . This value of the color excess 
was found from the intrinsic [U — B] color and reddening- 
free Geneva X and Y parameters (Cramer! 1 1982b . This esti- 
mate agrees with the E(B - V) r ange of 0.04- 0.08 that one 
can infer from the dust maps of iLuckd d 19781) and the re- 
vised Hipparcos par allax of HD 137509, n = 5.12 + 0.38 mas 
(I van Leeuwenl2007l) . The observed Stromgren and Geneva pho- 
tometric parameters of HD 137509 w ere t aken from the cata- 
logues of lHauck & Mermilliodl d 1998b and lRufenerj (1 19891) . re- 
spect ively. The observed value of the Aa parameter was adopted 
from lPaunzen et all d2005l) . These photometric characteristics of 
HD 137509 are summarized in the first row of Table [2] which 
also shows the results of the synthetic color calculations. We 
present results only for those models which fit the HB and Hy 
line profiles reasonably well. We have also looked into possibil- 
ity of using infra-re d photometric measureme nts (e.g., 2MASS0 
or observations by Groote & Kaufmann] d!983l) ) to constrain the 
model parameters. However, the difference between all models 
fluxes of HD 137509 considered below is less than the observa- 
tional uncertainty in this spectral range. 

There are several things to note in Table|2] First, it is evident 
that the change in photometric parameters associated with the 
introduction of the magnetic field and individual abundances in 
the model atmosphere calculations amounts to a few hundreds of 
a magnitude compared with the scaled-solar abundances model 
and that, generally, magnetic models allow us to improve the fit 
to the observed values of all photometric indices. The only ex- 
ceptions are the parameter Y for the model with T e g = 12 250 K, 
logg = 3.8 (its value for non-magnetic model is closer to the 
observed one compared to the magnetic model) and the index 
B2 - B (which tends to increase if the magnetic field is intro- 
duced to calculations, while decreasing is needed to fit the ob- 
servations). Although the overall changes in the color indices are 
not very large even for such a strong magnetic field, we conclude 
that taking into account Zeeman splitting and magnetic intensifi- 
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cation of spectral lines clearly is a necessary ingredient for mod- 
elling photometric observables of magnetic stars. Of course, the 
dependence of photometric colors on the magnetic field is sensi- 
tive to the chemical composition of the star. The increased metal 
content of the stellar atmosphere has a large impact on the photo- 
metric colors and thus HD 137509 represents one of the extreme 
cases due to the strong overabundance of iron-peak elements in 
its atmosphere. 

To distinguish between the effects of magnetic field and indi- 
vidual abundance we also compare synthetic colors for the mod- 
els calculated with the same fundamental parameters (T e ff = 
12750K, logg = 3.8) but different assumptions about abun- 
dances and magnetic field. Comparing the models with scaled- 
solar abundances with individual abundances model and the 
latter model with the magnetic model atmosphere+individual 
abundances, one can infer that the effect of the very strong mag- 
netic field on photometric observables is generally comparable 
in magnitude with the effect of individual abundances alone. 
Interestingly, for some of the color indices the difference be- 
tween the models calculated with and without magnetic field 
is larger than the same difference between the non-magnetic 
model with individual abundances and the reference scaled-solar 
model. This is the case for c\ , Aa, U -B,X,Y and Z. One should 
not forget that we are dealing with an effect depending on many 
parameters, i.e. the overall picture will depend upon the effective 
temperature and the abundance pattern, and thus could differ a 
lot from one star to another even if the field strength is the same. 
Limiting our conclusions to HD 137509, we find that the effect 
of the magnetic field is of no less importance than the effect of 
individual abundances for the majority of the photometric colors 
considered in the present study. 

Finally, we have also investigated the importance of accurate 
photoionization cross-sections for individu al states of Si 1 an d 
Sin ions taken from the TOPBase database dCunto et al.ll 1993b - 
The importance of accurate photoionization cross-sections wit h 
full resonance structure of Si 11 was noted bv lLanz et alj (Q996). 
They found that accounting for wide TOPBase resonances in the 
spectrum synthesis calculations leads to a good fit to the ob- 
served UV spectra of the Bp star HD 34452 (T eS = 13 650 K, 
\ogg = 4.0), which has an enhanced Si abundance. This could 
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Table 2. Observed and reddening-corrected calculated photometric parameters of HD 137509. Theoretical values were calculated 
for the models with ("mag") and without ("non") magnetic field, using individual abundances. Predictions for the scaled-solar 
abundance non-magnetic model of Paper I are given for comparison. The last row reports predictions for the magnetic and non- 
magnetic models with revised abundances determined in the present paper (see text) 





b-y 


m\ 




Aa 


U-B 


V-B 


G-B 


B { -B 


B 2 -B 


Vi~B 


X 


Y 


Z 


observations 


-0.095 


0.183 


0.411 


0.066 


0.761 


1.136 


2.294 


0.811 


1.566 


1.836 


0.762 


0.076 


-0.067 


tl2250g38, mag 


-0.054 


0.142 


0.640 


0.057 


0.995 


1.088 


2.265 


0.838 


1.574 


1.801 


1.026 


0.122 


-0.045 


tl2250g38, non 


-0.040 


0.134 


0.653 


0.033 


1.010 


1.069 


2.261 


0.846 


1.564 


1.777 


1.046 


0.073 


-0.028 


tl2750g40, mag 


-0.057 


0.141 


0.572 


0.054 


0.923 


1.094 


2.277 


0.838 


1.574 


1.806 


0.937 


0.088 


-0.044 


tl2750g40, non 


-0.044 


0.135 


0.586 


0.031 


0.941 


1.075 


2.271 


0.846 


1.565 


1.783 


0.957 


0.045 


-0.029 


tl3250g40, mag 


-0.059 


0.135 


0.516 


0.051 


0.860 


1.105 


2.293 


0.834 


1.578 


1.815 


0.863 


0.077 


-0.041 


tl3250g40, non 


-0.046 


0.130 


0.529 


0.030 


0.876 


1.087 


2.285 


0.841 


1.570 


1.794 


0.880 


0.040 


-0.028 


tl3750g42, mag 


-0.062 


0.138 


0.455 


0.046 


0.798 


1.106 


2.299 


0.839 


1.574 


1.815 


0.780 


0.029 


-0.040 


tl3750g42, non 


-0.051 


0.138 


0.468 


0.030 


0.815 


1.081 


2.282 


0.850 


1.561 


1.788 


0.790 


-0.023 


-0.030 


tl2750g38, mag 


-0.056 


0.135 


0.573 


0.054 


0.919 


1.101 


2.283 


0.833 


1.579 


1.812 


0.936 


0.107 


-0.043 


tl2750g38, non 


-0.043 


0.129 


0.588 


0.032 


0.936 


1.082 


2.277 


0.841 


1.570 


1.790 


0.957 


0.066 


-0.028 


tl2750g38, scaled 


-0.026 


0.116 


0.581 


0.014 


0.932 


1.061 


2.264 


0.847 


1.564 


1.766 


0.942 


0.032 


-0.016 


*tl3750g42, mag 


-0.081 


0.153 


0.453 


0.062 


0.800 


1.136 


2.321 


0.829 


1.585 


1.847 


0.799 


0.081 


-0.051 


*tl3750g42, non 


-0.064 


0.144 


0.470 


0.037 


0.820 


1.111 


2.307 


0.837 


1.576 


1.819 


0.818 


0.040 


-0.036 



be an important effect for HD 137509 because this star is also 
characterized by an excess of Si in its atmosphere. Thus, we 
have implemented available Si i and Si n cross-section data from 
the TOPBase web pagf] in the LLmodels atmospheric calcula- 
tions. In this modelling we did not use any additional techniques, 
such as resonance-averaging of photoionization cross-sections 
data or combining energy levels to superlevels which are usu- 
ally suggested to reduce computational expenses. All the en- 
ergy levels for both ions were explicitly included in the opacity 
calculations. We found that even for the hottest model consid- 
ered here the effect of detailed photoionization cross-sections is 
negligible for the majority of color indices in comparison with 
the bound-bound treatment of resonance structure via autoion- 
izing lines. The strongest deviation between the model with the 
TOPBase data and the model with Si i and Si n co ntinuous opac- 
ities taken from the Atlas 12 code dKuruczlll993bl) is found for 
ci (-0.004mag), U-B (-0.007 mag), X (-0.01 mag) and Y 
(-0.01 mag). We found that the main reason for such a small 
influence of the detailed photoionization data is that some of the 
strong autoionizing lines of Sii and Sin are already included 
in the VALD database and hence in the standard LLmodels input 
line list. This means that the photoionization cross-sections from 
Atlas 12 together with the autoionizing lines from VALD pro- 
vide a reasonably good approximation to the opacity calculated 
from the more accurate TOPBase data, at least for temperatures 
around 12 750 K. It is worth to note that here we are mainly inter- 
ested in the color indices, however, the accurate photoionization 
cross-sections are essential when it comes to the comparison of 
theoretical calculations with high-resolution UV spectra of Si- 
rich stars. 

3.5. An improved model for HD 137509 

The behavior of b - y, c\ and X indices suggests that the effec- 
tive temperature of HD 137509 could be higher than the initially 
derived value of 12 750 K which was used in Paper I. To im- 
prove the model parameters in a consistent way we redetermined 
abundances of chemical elements using model atmosphere with 
T e s = 13750K, logg = 4.2 calculated with magnetic field and 
individual abundances from the second column of Table Q] The 
resulting new abundances, determined using the same method- 



4 http : //vizier . u-strasbg . fr/topbase/topbase . html 



ology and spectral regions as in Paper I, are listed in the third 
column of Table [TJ The new chemical abundances are systemat- 
ically higher by 0.2-0.4 dex compared to the values reported in 
Paper I. This is mostly a result of using the model atmosphere 
with 1000 K higher effective temperature. The atomic lines used 
for abundances analysis become systematically shallower with 
the increase of atmospheric temperature thus resulting in pos- 
itive abundance corrections. Using these new abundances, we 
made a next iteration by re-calculating model atmosphere to 
ensure consistency of the model temperature-pressure structure 
and adopted abundances. This final model allowed us to fit the 
hydrogen line profiles as good as previous models. It also im- 
proves, with a few exceptions, the agreement between observed 
and computed photometric parameters as demonstrated by the 
two last rows of Table [2] (models marked with asterisks). Thus, 
we can conclude that the effective temperature of HD 137509 
should be around 13750K. 

We would like to note that exact fitting of the photometric 
indices is hardly feasible for any model due to additional un- 
known parameters involved in the modeling process. One of the 
effects which we have ignored is inhomogeneo us surface dis- 
tributi on of chemical elements in HD 137509. Mathvs & Lanz 
(119971) discussed rotational modulation of the equivalent widths 
of Cr ii, Si ii and Fe n lines, and the light variability in the U, B 
and V filters with the maximum amplitude of about 0.05 mag for 
the U filter. These phenomena are usually attributed to the pres- 
ence of abundance spots on the stellar surface. Therefore, to fully 
characterize atmospheric properties of HD 137509, one would 
need to obtain time-resolved observations and reconstruct sur- 
face distribution of magnetic field and c hemical abundances with 
the Doppler imaging technique (e.g.. iKochukhov et alj 12004). 
Furth ermore, the vertical str atification of chemical elements 
(e.g., Kochuk hov et al.l 12006) could also take place in the at- 
mosphere of HD 137509. A refined analysis of these chemical 
inhomogeneity effects is outside the scope of our study. 

4. Summary and conclusions 

In the present paper we have constructed advanced theoretical 
stellar model atmospheres incorporating accurate treatment of 
the individual abundances pattern, Zeeman splitting, polarized 
radiative transfer and compared the results with the observations 
of extreme magnetic CP star HD 137509. With the mean surface 
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field of (B) » 29 kG, this object has the second largest magnetic 
field among CP stars. Strong overabundance of iron-peak ele- 
ments and extreme underabundance of helium in the atmosphere 
of this star opens a possibility to investigate the importance 
of taking individual abundances into account when construct- 
ing model atmospheres of magnetic chemically peculiar stars. 
Theoretical model atmosphere calculations were compared with 
the hydrogen line profiles and metallic spectrum of HD 137509. 
The Stromgren and Geneva photometric parameters were also 
investigated. The main conclusions of our study can be summa- 
rized as follows: 

- We found that the effect of individual abundances dominates 
the change of Balmer H/3 and Hy line profiles compared to 
the model with solar-scaled abundances. This implies that 
once the abundances of the star have been determined us- 
ing an approximate model, it is necessary to recalculate the 
model atmosphere in order to ensure the consistency be- 
tween abundance pattern and the model structure. The mag- 
netic field has less influence on the hydrogen lines, however 
it should be taken into account for stars with very strong 
magnetic fields. 

- Modification of the atmospheric temperature-pressure struc- 
ture due to presence of the magnetic field and peculiar abun- 
dances generally has little impact of the metal line profiles 
comparing to that of non-magnetic one. 

- For HD 137509, the effect of magnetic field on photometric 
colors is very important for some photometric parameters. 
Occasionally the combined impact of the magnetic field and 
the realistic chemistry is more important than the effect of us- 
ing only individual abundances. Generally, magnetic model 
atmospheres allowed us to obtain a better agreement between 
almost all observed and theoretical color indices. Thus, we 
can conclude that magnetic field should be taken into ac- 
count in the analysis of stars with strong magnetic fields. 

- We showed that the analysis of the spectra of such extreme 
Bp stars with strong magnetic fields and unusual chemistry 
as HD 137509 requires a self-consistent approach. Once the 
abundances of the most important elements are derived us- 
ing an approximate model atmosphere, it is necessary to re- 
compute the model with new abundances trying to fit various 
observables such as hydrogen line profiles, photometric col- 
ors, and energy distribution (if available). For HD 137509 we 
found that the simultaneous fit to the hydrogen line profiles 
and photometrical indices employing the model with both 
magnetic field and individual abundances included requires 
as much as 1000K correction to the effective temperature 
and 0.4 dex correction to the surface gravity compared to the 
results obtained using simple scaled-solar models. 

- We expect that the overall energy distribution of the star 
is strongly modified by magnetic line blanketing. However, 
the lack of accurate energy distribution covering meaningful 
wavelength region in the spectrum of HD 137509 precludes 
us from reaching robust quantitative conclusions. We em- 
phasize that availability of the flux distributions for this and 
other magnetic CP stars is of great importance for the stel- 
lar parameter determination and for verification of the new 
generation model atmospheres of magnetic CP stars. 
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